!***************************************
FINISH
/CLEAR
/FILN, APDL_EXAMPLE  
/TITLE, VORTEX
/UNITS,SI                 !ùʵλ

!*****************************************
PI=3.14159265
RR=8.75E-3                !߻Բľ
RB=4.2176E-3              !߻Բ뾶
ALAF=30                   !߷
HW=0.03                   !ݸ
WO=0.0835                 !̰뾶
H=0.02                    !̺
EXX=2.07E11               !òϵĵģ
NUXYQQ=0.3                !òϵĲɱ
GXYQQ=0.5*EXX/(1+NUXYQQ)  !òϵļеģ
DENSQQ=7850               !òϵܶ
ALPQQ=9.5E-6              !òϵϵ
LAYER1=6                  !Ԫ

!****************ǰ***************************
/PREP7  
!************趨Ԫͼ***************
ET,1,95                   !1ŵԪ-SOLID95
MP,EX,1,EXX               !趨ϵĵģ
MP,GXY,1,GXYQQ            !趨ϵļеģ
MP,NUXY,1,NUXYQQ          !趨ϵĲɱ
MP,DENS,1,DENSQQ          !趨ϵܶ
MP,ALPX,1,ALPQQ           !趨Xϵ
MP,ALPY,1,ALPQQ           !趨Yϵ
MP,ALPZ,1,ALPQQ           !趨Zϵ


*AFUN,DEG                 !ǺöǵĴС

!ͨؼ
!ȡȦڲཥߵĹؼ
*DIM,INNER,ARRAY,111,2                                 !111X2Ķά
*DO,II,0,1100,10                                       !趨ѭII0仯1100Ϊ10
INNER(II/10+1,1)=RB*(COS(II)+(II+ALAF)*PI/180*SIN(II)) !ȦڲཥߵĹؼX
INNER(II/10+1,2)=RB*(SIN(II)-(II+ALAF)*PI/180*COS(II)) !ȦڲཥߵĹؼY
*ENDDO                                                 !ѭ
!*STATUS,INNER,1,111,1,2,1                             !бʾ

!ȡȦཥߵĹؼ
*DIM,OUTER,ARRAY,108,2                                 !108X2Ķά
!ע⣬ȦıںȦڡཥʼƫƫΪ30ȣII30ʼ仯
*DO,II,30,1100,10                                      !趨ѭII30仯1100Ϊ10
OUTER(II/10-2,1)=RB*(COS(II)+(II-ALAF)*PI/180*SIN(II)) !ȦཥߵĹؼX
OUTER(II/10-2,2)=RB*(SIN(II)-(II-ALAF)*PI/180*COS(II)) !ȦཥߵĹؼY
*ENDDO                                                 !ѭ


!**********************ģ************************
!****************֪Ŷؼ********************
*DO,II,0,1100,10                                       !趨ѭII0仯1100Ϊ10
K,II/10+1,INNER(II/10+1,1),INNER(II/10+1,2)            !ڲཥ߹ؼ
*IF,II,GE,30,THEN                                      !жϵѭII=30֮
K,II/10+201,OUTER(II/10-2,1),OUTER(II/10-2,2)          !ཥ߹ؼ
*ENDIF                                                 !жϽ
*ENDDO                                                 !ѭ


!************⿪߲䶨3㣬Աߵķ***********
*DO,FAI,30,50,10/3
X2=RB*(COS(FAI)+(FAI-ALAF)*PI/180*SIN(FAI))         !߷õX
Y2=RB*(SIN(FAI)-(FAI-ALAF)*PI/180*COS(FAI))         !߷õY
K,FAI*3/10+192,X2,Y2                                !ؼ
*ENDDO                                              

!*************ɵ㹹*********************************
*DO,AAK,1,106,5                                     !趨ѭΪ5ʼΪ1ֹ106
BSPLIN,AAK,AAK+1,AAK+2,AAK+3,AAK+4,AAK+5            !ڲཥ
*ENDDO                                              !ѭ

*DO,AAL,201,306,5                                   !趨ѭΪ5ʼΪ201ֹ306
BSPLIN,AAL,AAL+1,AAL+2,AAL+3,AAL+4,AAL+5            !ཥ
*ENDDO                                              !ѭ

!************ߣΪ׼*******************
*DO,JJ,1,111,5                                      !趨ѭΪ5ʼΪ1ֹ111
L,JJ,JJ+200                                         !ߣΪ׼
*ENDDO                                              !ѭ



!************̽ߵĩɵ**************
LFILLT,22,67,0.001
LFILLT,67,44,0.001
!*************ɵһ߶*******
LCOMB,68,67                       !67ԭ߶45Ϊһ߶                     
LCOMB,69,67                       !46ԭ߶45Ϊһ߶   
L2TAN,-1,23                       !߶233еֱ
LDEL,45

!**********************߹****************
*DO,ABC,2,22,1 
AL,ABC,ABC+45,ABC+22,ABC+44
*ENDDO
AL,23,68,1,46

!**********ɾĹؼ**********************
ALLSEL                      !ȫѡ
LSLA                        !ѡϵ
KSLL                        !ѡϵĵ
KSEL,INVERT                 !ѡ
KDELE,ALL                   !ɾ

ALLSEL
*GET,MAXNUM,KP,,NUM,MAX     !ؼ

!*************ɵZ****************
K,MAXNUM+1                    !ؼ
K,MAXNUM+2,,,HW               !ؼ㣨ϵXYһ£ZΪHW
L,MAXNUM+1,MAXNUM+2           !ͨ㶨߶
LESIZE,45,,,LAYER1,2          !¶߶趨ַ 
VDRAG,ALL,,,,,,45             !¶ֱ

NUMMRG,ALL                    !Ա
NUMCMP,ALL                    !ѹ


CYLIND,0,WO,,-H,,360          !
ALLSEL                        !ȫѡ
VGLUE,ALL                     !е

!*******************************
TYPE,1                        !趨Ԫ1
MAT,1                         !ָ1
ESIZE,0.006,0,                !趨ֵԪС
VSEL,ALL                      !ѡе
VSEL,U,VOLU,,23               !ȥ
VSWEEP,ALL                    !ӳ񻮷

MSHAPE,1,3D 
MSHKEY,0
VSEL,S,LOC,Z,0,-H               !ȥ
VMESH,ALL